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en Abstract 

^D We introduce an algorithm for the pricing of finite expiry American options 

driven by Levy processes. The idea is to tweak Carr's 'Canadisation' method, cf. 
O I Carr [9] (see also Bouchard et al [5]), in such a way that the adjusted algorithm is 

^^^ viable for any Levy process whose law at an independent, exponentially distributed 

\^ time consists of a (possibly infinite) mixture of exponentials. This includes Brow- 

^~~^ nian motion plus (hyper) exponential jumps, but also the recently introduced rich 

class of so-called meromorphic Levy processes, cf. Kyprianou et al [1,6] . This class 
contains all Levy processes whose Levy measure is an infinite mixture of exponen- 
tials which can generate both finite and infinite jump activity. Levy processes well 
known in mathematical finance can in a straightforward way be obtained as a limit 
of meromorphic Levy processes. We work out the algorithm in detail for the classic 
example of the American put, and we illustrate the results with some numerics. 
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^ 1 Introduction 

j> Let X be a Levy process defined on a filtered probability space (fi, J-", F,P), where F = 

k> {J^t)t>o is the filtration generated by X which is naturally enlarged (cL Definition 1.3.38 

in |1]). Recall that a Levy process is characterised by stationary, independent increments 
and paths which are right continuous and have left limits, and its law is characterised by 
the characteristic exponent "^ defined through E[e'^'''^'] = e~^^^^^ for all t > and z G M. 
According to the Levy-Khintchine formula there exist a > 0, a G M and a measure 11 (the 
Levy measure) concentrated on M \ {0} satisfying f^{l A x^) Il{dx) < oo such that 
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^i{z) = ^z^ + iaz+ / (1 - e''"" + l{|^.|<i}izx) U{dx) 

for all 2; G M. The tuple (a, a, 11) is usually refered to as the Levy triplet. For a further 
introduction see e.g. the textbooks [3| [T8l l2l] . 

For a; G M we denote by P^. the law of X when it is started at x and we abbreviate 
P = Pq. Accordingly we shall write E^. and E for the associated expectation operators. 
We denote by T the set of F-stopping times. 

Consider the following classic optimal stopping problem with expiry date T G [0, 00): 

v{T,x) = supE, [e-^(^^^)/(X.^T)] , (l-l) 

where / is a payoff function and r > the discount rate. Such problems have been studied 
thoroughly and have applications in many fields, where option pricing theory in modern 
mathematical finance is maybe one of the most appealing ones. See e.g. the textbook [23j 
for a recent overview of the general theory concerning optimal stopping problems, and 
e.g. [I2] for an introduction to mathematical finance and explanation of the related terms 
used here. For example, if in a typical Black & Scholes setup the price of the risky asset 
5* is assumed to evolve as St = exp(Xt) for all t > while the price of the riskless asset 
S^ is given by S^ = e*"* for all t > then v can be interpreted*] as the fair price of an 
American option on the risky asset which expires after T time units and x = Xq = log 5*0. 
One of the classic examples is the American put option, which is characterised by the 
payoff function f{x) = {K — e^)^ for some K > 0. Here a^ := max{a, 0}. Such Levy 
process driven market models have received considerable attention in recent years, in an 
attempt to overcome some of the limitations of the classic Black & Scholes model where 
X is a Brownian motion. See e.g. [H] and the references therein. 



It is well known that in general no closed form formula exists for (1.1). Results do 
exist for certain combinations of payoff functions / and Levy processes X when T = 00, 
but when T < 00 — as we assume throughout this paper — one typically has to settle 
for approximating v rather. 

In his celebrated paper |9] (see also [5j), Carr proposes a method for this he refers to 
as 'Canadisation'. The idea is to introduce a 'stochastic time grid', which is a grid where 
the distances between consecutive grid points forms a sequence of independent random 
variables with an identical exponential distribution. The grid gets finer as the common 
mean of the exponential distributions vanishes. The value of an American option with 
expiry date T is then approximated by the value of the same American option but with 
the expiry date replaced by a grid point close to T. This approximating value function 
can be determined by performing a backwards induction over the grid points. Carr works 



*In addition it is required that S/S'^ is a martingale, li X is not a Brownian motion there are typically 
multiple measures Q equivalent to each other under which S/S^ is a martingale. For this paper it is 
not relevant which measure is used, as long as X is a Levy process under the chosen measure and P is 
understood to be the chosen measure 



out this algorithm for the American put option driven by a Brownian motion. See also [8] 
and the references therein for more general applications of Carr's algorithm. 

As shown in [§| Carr's Canadisation is also viable for Brownian motion with exponen- 
tial jumps. However, it does not seem easy to apply this algorithm to more general Levy 
processes due to the fact that the backwards induction over the grid points is involved 
and (hence) imposes restrictions. This paper introduces a variation of the Canadisation 
algorithm, where some changes to the original setup are made so that the backwards in- 
duction over the grid points becomes more straightforward. Indeed, the resulting adjusted 
algorithm allows to approximate v by elementary functions not just for Brownian motion 
plus exponential jumps but for any so-called meromorphic Levy process. 

The class of meromorphic Levy processes was recently introduced in [16] and can be 
defined as the family of all Levy processes whose Levy measure has a density that consists 
of an infinite mixture of exponentials, see Section |4]for further details and examples of such 
processes. This is a rich class, where in particular also infinite activity jumps are possible. 
It contains for instance the family of so-called /3-processes (cf. [E]). Furthermore many 
Levy processes that play a prominent role in mathematical finance can in a straightforward 
way be obtained by considering limiting cases of meromorphic Levy processes, see |T3] and 
[IB] . Examples include generahzed tempered stable processes (cf. [H]), Kobol processes 
(cf. j6] and [7]), tempered stable processes (cf. [2]) and CGMY processes (cf. p^). 

In particular meromorphic Levy processes have the property that the law of the pro- 
cess evaluated at an independent, exponentially distributed time is an infinite mixture 



of exponentials (cf. Proposition 4.2). (In fact meromorphic Levy processes are the only 
Levy processes for which this is the case.) Note that for a Brownian motion with hyper- 
exponential jumps this law is a finite mixture of exponentials. This property is key in our 
algorithm, as becomes clear in Section |4j Indeed, if this is the case every step in the algo- 
rithm can be worked out explicitly in terms of elementary functions provided the payoff 



function / in (LI) can (piecewise) be expressed as a linear combination of functions of 
the form Ax'^e^^ + C for A,B,C G M and i G N. This is for instance the case for the 
classic example of the American put option, where f{x) = {K — e^)^ for some K > 0. For 
any payoff function not of this type it is a straightforward exercise to approximate it by 



a function of this type and estimate the error due to this approximation, cf. Remark 4.5 
Pricing American options driven by Levy processes has received considerable attention 
in recent years. Indeed several other possible techniques were developed, including (but 
not limited to) tree-based methods (see e.g. [21]) and variational methods (see e.g. [22]). 
See also the overview in [20] . A detailed comparison of all available methods with the one 
proposed in this paper would constitute a paper in itself and we do not embark on such a 
task here. The authors feel the main appeal of the method presented in this paper lies in 
the fact that it is a 'light weight' method which is probabilistic in nature and does not need 
any involved proofs or tools from other fields. Furthermore the computer implementation 
is rather straightforward as only elementary operations are required (together with the 
straightforward operation of numerical root finding). Finally the method is fiexible and can 



easily be adjusted to deal with situations where for instance the option has a Bermudan 
type structure, i.e. where exercising is only allowed in certain subsets of [0,T], or there 
is some path dependency in the payoff (for instance when the payoff also depends on the 



path supremum). Cf. Remark 4.5 



The rest of this paper is organised as follows. In Section [2] we introduce the algorithm 
in detail and spend some time discussing the difference with the original Canadisation 
approach. In Section [3] we make the algorithm rigorous. In Section |4] we work out the 
algorithm in detail for the prominent example of the American put, where the driving 
process is a meromorphic Levy process (or a Brownian motion plus (hyper)exponential 
jumps). In that section we also provide some more details concerning meromorphic Levy 
processes. We find that the functions approximating v generated by our algorithm have 
easy to implement, explicit formulae. We conclude by presenting some numerical results 
in Section m 

2 The algorithm and comparison with Canadisation 

Let us introduce the algorithm in detail. For any n G N, enlarge the above probability 
space to contain a Poisson process A^(") with intensity n, independent of X. Denote the 
fc-th jump time of A^(") by T^ , i.e. 

T^") := inf{t > I Ar,(") > k} (2.2) 

with Tq := 0. Let F*^") denote the naturally enlarged filtration generated by the pair 
(X, A^*^"-*). Note that for any T > 0, if {k{n))n>i is a sequence of natural numbers such 
that k{n)/n — )■ T as n — )■ cxd then TyK — )■ T a.s. by the law of large numbers. 

Furthermore, denote by T*-"-* the set of F'^"^-stopping times on this enlarged probability 
space that only take values on the grid points {0 = Tq , T|" , . . .}, that is 

f(") := |r r is an F(")-stopping time & r G {0 = T^''\t[''\ . . .} | . 



Now consider modifying the original optimal stopping problem (LI) in three steps. 
First, replace the deterministic expiry date T by the random variable T^ (for suitably 
chosen k). Second, replace the set of stopping times T over which is optimized by T*-"-*. 
Finally, replace the discount factor e"*""^ by D^"\t) defined as 

oo 

D(")(r):=5^1^^^^(.,^e-"/". (2.3) 

j=0 

Together this amounts to defining for each ra G N and A; G N: 

Vt\x) := sup E,. bW(r A T^ )/K,^("))l • (2-4) 



The usefulness of this setup rehes on the following two facts. For any expiry date 
T > 0, choosing a sequence (fc(n))n>i such that k{n)/n - ^ T as n — )■ oo we have for any 
X G M that V^'^\{x) — !■ v{T,x) as n — )■ oo (cf. Theorem 

n G N the sequence of functions (V^"' )k>o satisfies the foL 



3.1 



(ii)). Furthermore, for any 



owing recursion: 



Ko(")(x) 



fix) 



Vt\x) 



max 



/(x),e->E,. C1(X^(.)) 



for /c > 1, 



(2.5) 



where ^^'^' is an exponentially distributed random variable with mean 1/n, independent 



of X. This is easily seen by a dynamic programming argument. Indeed, consider (2.4) for 



some X G M. The available stopping times allow to stop either immediately or at T{ or 
later. Stopping immediately yields a payoff of f{x) while the (current) value of waiting 
until T}" or later equals the discounted expected value of the option with expiry date 



T^ -^ for X 



X^(n). Hence Vf^ 



(«), 



x) equals the maximum of the two. 



The main point of this setup is that the recursion (2.5) is very straightforward. Indeed 



as already alluded to in the Introduction, for any meromorphic Levy process the law of 



X^(n) is a mixture of exponentials. As a consequence (2.5) can be worked out explicitly 



for example if f{x) = {K — e^)~^, where first the expectation is computed explicitly and 
next the maximum is easily determined, cf. Section |4j 

Now let us for comparison recall the original Canadisation algorithm introduced by 



<n) 



Carr. With n ' as defined in (2.2), define for all n G N and fc G N 



r(n) 



V;^"^(x) = supE,. 



-r(rAT(")) 



fix^ 



AT 



("), 



where the supremum is taken over all F'^'^^-stopping times. It can be shown (cf. also ^5]) 
that this yields the following recursion: 



r{n), 



-rg(")f>(n) 



V^'ix) = supE. l^,<^,„)}e-^V(Xr) + l{.>5(n)}e-'-«^ Vn(X5w) 



SUpEa; 
reT 



-("+")V(Xr) + n I e" ("+")" V;^"J(X„) du 



(2.6) 



To work out this recursion, general theory of optimal stopping (cf. e.g. [23]) dictates that 



the optimal stopping time in (2.6) is the first time that X enters a subset S of the state 
space M, say r^. Hence, in order to find an explicit expression for (2.6) one requirement is 
to know the mean of e~^^~^"'^'^^ /{Xt-^). Carr worked with Brownian motion in which case 
this quantity is available in closed form. However once X has jumps this becomes much 
more involved. For meromorphic Levy processes expressions are available that could in 
principle be used when iS is a halfline or an interval (cf. Theorem 3 and Theorem 5 in |16j ) 
but these are quite involved and complicated to implement. If S consists of a union of 
intervals the situation becomes even much more complicated. Namely, for the overshoot 
the possibility that X jumps between different disconnected parts of S'^ (the complement 



of S) would have to be taken into account. That is, for some x G iS^ the rhs of (2.6) would 



depend on the values the function V^ takes in intervals that are part of S'^ and that X 
could reach by a jump without entering S first. 

Furthermore, in the case of an American put, i.e. f{x) = {K — e^)^, it can be shown 
that S = (—00, a;*] for some x* (depending on n and k) and it remains to determine 
X* . For other payoff functions it may a priori very well not be clear at all what shape S 
exactly takes, which makes the above task even more challenging. 



On the other hand, the recursion (2.5) generated by our algorithm is much more 



straightforward to implement. In the case of the American put it can a priori be shown 



that a level x* (depending on n and k) exists such that (cf. Proposition 4.3) 



vt\^) 




V^fc-i(^€("^ 



a X < X* 
a X > X* 



and working out this recursion is easily possible for any meromorphic Levy process. For 



other payoff functions where a priori it is not clear how the maximum in (2.5) works out 



it is still rather straightforward to let a computer do the work. The first step would be to 
determine an explicit formula for 



X ^-)■ e' 



-'■/"E^ 



V^'^^{X^(n)] 



(2.7) 



.(n). 



and then V^ (x) is the maximum of (2.7) and f{x). With formulae available for both (2.7) 



and / this is rather straightforward to implement using a computer. See also Remark 4.5 



n . 



Remark 2.1. A natural question is whether it would not be better to use a classic de- 
terministic grid rather than a stochastic grid. That is to say, setting tk = kT/n for 
k = 0,1, . . . ,n one could consider approximating v by narrowing the set of stopping times 
over which is optimised to the set of stopping times taking values in {to,ti, . . . ,t„} only. 
Denoting the resulting value function with expiry date tk by U]^ the sequence of functions 



(ut^: 



fc=o,i,...,n ^s determined by the recursion 



Ut\x) 



/(^) 



Ut\x) 



max |/(x), e-^/"E, [utMx^/^)] ] for k = 1 



,n. 



(2. 



The problem is that this recursion involves the law o/Xi/„. Even if this is known (which 
generally is not the case except for a few notable exceptions) it is typically not of a friendly 



enough nature to work out the recursion (2.8) explicitly. Indeed consider for instance a 
Brownian motion. 



3 Main result & proof 



This section is dedicated to making the algorithm outlined in the above Section[2]rigorous. 



Theore m 3. 1. Let f be a bounded and continuous function. Let the valu e fun ction v be 
given by (1.1). For each n > 1 and k > let the function V^ be given by (2.4) ■ We have 
the following. 

(i) For any n>l, the sequence of functions (V^ )fc>o satisfies the following recursion: 



Vt\x) = fix), Vt\x) = max |/(x), e-^/^E, [v^^^ll^^w)] } for k > 1, (3.1) 

where ^^"'^ is independent of X and follows an exponential distribution with mean 
1/n. 

(a) For any T G [0, oo), if {k{n))n>i is a sequence such that k{n)/n —)■ T as n —t- oo 
then: 



T, 



tin)-^T a.s. and V;^"J)(x 



v{T, x) for all x G M. 



Proof. Ad (i). This is just the dynamic programming principle applied to (2.4). Indeed 
for any k>l 



r{n) 



vr{^) 



sup Ea; 

rer(") 



L{-o} + l^,>^(.)j) D^''\r A Tf )/(^.,^.o)' 



sup Ea; 



sup E. l|.=o}/(x) + 1. >^(.).E,. D^-\t ^T^)f{X^^^,^,) 



l{r=,}f{x) + l^^>^(,0}e-^/'^E,. \vt\{X, 



'-^{n), 



max |/(x), e-^/"E, [v^^-ll^gw)] } • 



■Jrpin) 

1 J J 



Ad (ii). Ty^K — )■ T a.s. as n — >■ oo is just an application of the law of the large numbers. 

For proving that V^^?\{x) — )■ f (T, x) as n — )■ oo for any a; G M we have to do a bit more 
work. In this proof we will also need an auxiliary sequence of functions that is generated 
by our algorithm if we leave out the step of adjusting the discounting factor, that is, let 
for n > 1, A; > the function VJ: be given by (compare with (2.4)): 



Vt\x) 



sup Ej. 

rer(") 



n(") 



rirAT}"') f(v \ 



Take any a; G M. The result follows from the following three steps. 
Step 1: 



rin) 



Step 2: 



limsup V2(„)(x) < v{T, x). 



lim inf Kvll (x) > v{T, x) . 



(3.2) 

(3.3) 
(3.4) 



Step 3: 



lim 

n—^OD 



vS^i-)-V^ix] 



0. 



(3.5) 



Step 1. We appeal to Carr's Canadisation as outlined in Section [2} however we follow 
the setup from [5]. Define for n > 1 the sequence of functions (U^ )k>o as follows (recall 
that T is the set of stopping times with respect to the completed filtration generated by 
X): 



rW/ 



rW/ 



Ur{x) = f{x), Ur {x) = snpE^ 



Jo 



It is shown in [5] that 



UyK{x) -^ v{T,x) as n -)■ cx). 



for A; > 1. 
(3.6) 

(3.7) 



By a dynamic programming argument the sequence (V^ )k>o as defined in (3.2) satisfies 



the recursion 



r{n) 



in). 



V.^'^ix) = fix), V^ix) = max /(x),E, 



-r^f") T>(n) 



^fc-l(^e("' 



for A; > 1, 



where ^*^"^ is an exponentially distrib uted random variable with mean 1/n, independent 



of X. Comparing this recursion with (3.6), we see that V"^_{ < f/^_i imphes V^: < U^ , 



since in (3.6) the supremum is at least as large as taking either r = or r = o o. Since 
yW _ jjw _ J j^ follows that V^ < Uj^ for all k,n, which together with (3.7) indeed 



yields (3.3). 



Step 2. Let e > be arbitrary and let us show that 



liminf Kvl(a;) > v{T,x) - e. 



(3. 



For this, first note that since X is continuous over stopping times there exists a r' taking 
values in = to < ti < . . . < tAT = T such that 

E, [e-'-^7(X.O] >v{T,x)-e. (3.9) 

Define an approximating sequence of stopping times for r' by setting for any n > 1 

Now we claim that as n — )■ oo: 



E, 



e— '"7(x.w; 



E, 



e-'-^ f{X,. 



(3.10) 



and 



E, 






0. 



Using this we indeed arrive at (3.8): 



(3.11) 



77, — ^no ^ ' 



lim inf E, 



lim inf E^ 



-K.Wat';'))) ..^ . 

fc{n) 



e— '"V(X,(.); 



E, 



e-^^ /(X. 



t;(T,x) -£, 



so let us prove (3.10) & (3.11). 



Ad (3.10): By construction of r' and a^"^"* we have 



E, 



e-^^7(^r' 



-E,. 



e— '"V(^.(.o; 



N 
j=0 



where 



af):=inf{Tf |Ti")>t,}, 
so it is enough to show that for each z = 0, . . . , A^ 



(3.12) 



J") 



e ■"« 7(^^(")) — ^ e "^'Vl^tJ as n -> oo. 
Since / is continuous and bounded, say by M, we have 



(3.13) 



J") 



e-^''^ /(XM)-e"^'7(^t. 



< 



J") 



e-^'^^ /(X^(„))-e-'^VK(")) 

rti tl V \ „-rti 



+ 



e-'^*7K("))-e-'-'7(^t. 

i 
(n) 



< M ( e-'"^ ' - e-"*' ) + e-'''' 



)+e-^'*'|/K("))-/(^tJ 



so that by uniform integrabihty and the continuous mapping theorem we have that (3.13) 
follows from 



a, — > ti and X («) 



X+- as n — )■ oo. 



(3.14) 



Now, recalling (|3.12|) and that Tc^ is a sum of i.i.d. exponentials with parameter n, due 



An) 



to the lack of memory property al — tj is equal in distribution to an exponentially 



distributed random variable with mean 1/n, independent of X, say ^'^"•*. So the first part 



of (3.14) is obvious and the second part is a direct consequence of stochastic continuity 
of X. 



Ad (3.11): as above, it suffices to show 



a 



(n) 



(a(")AT,^;:;))^0 and X^^ - X^(„,^^,„) ^ as n ^ oo. 



k(n) 



(3.15) 



For the first one, for any 5 > we have 



= F,. (a(") > T^;:;) + <5 & Tl-1^<T-6I2 

< Px (r,("i) <T- 5/2) + P, (a(") > T + 5/2) . 

The first probabihty on the last line vanishes as n —t- 00 on account of TyK — > T a.s. To 
see that the second probability vanishes as well, note that by construction r' < T so that 
by definition of a'-"^ the difference a^""^ — T is bounded above by inf {T^" | T^ > T} — T, 
but again by the lack of memory property this random variable follows an exponential 
distribution with mean 1/n. 



Finally, the second statement in (3.15) follows as above from the first statement and 
stochastic continuity of X. 

Step 3. Some standard manipulations, using Holders inequality and the boundedness 



of /, show that the definitions (3.2) and (2.4) yield: 



KS"1(x)-V,("\(x1 



k{n)^ 



k{n) ^ 



k(n) 

< M sup ^ E^ [1 

k{n) 

< M sup ^ E^. [ 



l^^^^^n)} 



l|^=r^")}J 



^-ri/n _ ^-rT. 



.(n) 



1/2 



E,^ 



~ri/n _ p-TT> 



(n) 



ni/2 
(3..16) 



Recalling that T> is the sum of i i.i.d. exponentials with mean 1/n, it is straightforward 
to compute 



E, 



g-"/n _ e^''^i 



(n) 



n 



in + 2r) 



- 2e""/" 



n 



{n + r)' 



+ e 



-2ri/n 



As this expression is increasing in i the supremum in (3.16) above is attained by r = Tj^K. 



Plugging this in it is clear from the vanishing variance of T^K (or directly from the above 



display) that (3.16) indeed vanishes as n — )■ 00. 



D 



10 



Remark 3.2. The assumption in the above Theorem 3J_ that f is bounded and continu- 
ous keeps the proof compact. However the boundedness can be weakened to a more liberal 
integrability condition, and continuity can be weakened to (at least) a discrete set of dis- 
continuities — provided compound Poisson processes are excluded. Details are left to the 
reader. 

4 Example: the American put driven by a meromor- 
phic Levy process 

In this section we work out in detail our algorithm for an American put option driven by 
a meromorphic Levy process. This is probably the most classic example of an option in a 
Black & Scholes type financial market. If the option has strike price K > and the price 
S of the risky asset on which the option is written is assumed to evolve as St = exp(Xf) 
for t > then the payoff function is given by 

f{x) = {K-eT- (4.1) 

Let us first briefly discuss meromorphic Levy processes in more detail, see also [16]. A 
function / : M>o — )■ M is said to be discrete completely monotone if 

f{x) = J2 ame"^"^ for all a; > 0, (4.2) 

m>l 

where CLm > 0, bm> for all m > 1 and 6^ — ^ oo as ?Ti — > oo. The definition is as follows. 

Definition 4.1. A Levy process X is meromorphic if the functions x ^-^ U (x) = 
n((x, oo)) and x ^-^ U (x) = n((— oo, —x)) are both discrete completely monotone. 

Alternatively, this is any Levy process whose Levy measure consists of an infinite 
mixture of exponentials on both the negative and positive halfline (cf. equation (5) in 
|16j). As mentioned in the Introduction, the key property of a meromorphic Levy process 
X that makes our algorithm yield explicit expressions is the fact that the law of X 
evaluated at an independent exponentially distributed random variable ^^''^ with mean 
1/q is a mixture of exponentials: 

/ N+ 7V_ \ 

P (X^(„ edx)=[ l|.<o} Yl cfU)e^-'^'^'' + l{.>o} E cL'Hj>^-'^^'^^' dx (4.3) 

V i=o i=o / 

for Ar_,iV^ e N U {oo} and sequences (cL'^(j)), (c?(j)), (C-^(j)), (C+Hj))> where for all 
(relevant) j: c^l\j) > 0, c^^\j) > 0, d'\j) < 0, ci\j) > 0; and 






i>o 



11 



Note that iV_ = iV_|_ = corresponds to Brownian motion (with drift), while iV_, iV+ G N 
corresponds to a Brownian motion with drift plus (one- or two-sided) so-called hyper- 
exponential jumps (a generalisation of exponential jumps). See e.g. [13] and [T]. For N_ = 
A^_i_ = oo this expression for the density holds if and only if X is a meromorphic Levy 
process (the proof relies heavily on the results in p5j). 

Proposition 4.2. Let X be a Levy process which is not compound Poisson. Then X 



is a meromorphic Levy process if and o nly i f X satisfies {4-3) for some (and then all) 



q > 0. Furthermore, if X indeed satisfies (4-3), then (C_ (j))j>o (resp. {(+ iJ))j>o) is the 



"(q), 



^(q), 



sequence of all real negative (resp. real positive) zeros of the function z \-^ q + '^(iz) (or, 
equivalently, of z \-^ ^{~^) ~ i)j and for all j > 



-4'(-d"(j))' ' n-c'-\j))' 

Here ^ is the characteristic exponent of X and $ the Laplace exponent. 

Proof. The 'only if part can be found in Kuznetsov et al. (TH], Theorem 2. Let us now 



assume that X satisfies (4.3) for some (and then all) g > 0. By using Theorem 1 in 



Kuznetsov et al. [IE], it is sufficient to show that the functions x i-> P(Xe(g) > x) and 



X I—)- P(Xe(g) < —x) are discrete completely monotone functions on ]R>o (recall (4.2)) 
Now 



oo jq_ JSf_ 



P(X^(,) >x)= ["£ c^l\3) eMC^\3)x) dx = 5^ -^^ exp(d'')(j» 



c 



(i) 



(^■) ^ a ^^A h. - _Ai)i 



If we now take aj = — j^j— > and bj = —(_ (j) > which is in fact increasing we see 

that indeed x H- F{Xe(g) > x) is a discrete completely monotone function. This can be 
proven for x i— )■ P(Xe(g) < —x) in the same way. D 

Let us briefly outline a few prominent examples of meromorphic Levy processes. The 
(3-class was introduced in [H]. It has a characteristic exponent given by 



2 

*(^) = ^z^ + iaz+^ {B{ai, 1 - Ai) - B{ai - iz/(3i, 1 - Ai)) 
z pi 

+ ^ (5(^2, 1 - A2) - B{a2 + iz/P2, 1 - A2)) 
P2 

where B{x,y) = r(x)r(y)/r(x + y) is the Beta function and the parameter ranges are 
a G M, o", Cj, ai, /3j > and Aj G (0, 3)\{1, 2}. The corresponding Levy measure has density 

^i-<oi''2^1_efcx)A2 + li->o}^i(i_e-/3i-)Ar 
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Note that setting a = and Aj G (0, 2) yields a process with paths of bounded variation, 
all other choices paths of unbounded variation. Furthermore infinite (resp. finite) activity 
in the jump component can be obtained by setting Aj G (1, 3) (resp. Aj ^ (1, 3)). 

A second example is the 0-class introduced in [15]. The Levy measure of such processes 
have a density of the form 

where G^ is the fc-th order (fractional) derivative of the theta function x ^-^• ^3(0,6"^): 
Qk{x) = ^^3(0, e-^) = l|fc=o} + 2 Y, r^2'=e-"'^ a; > 0. 

n>l 

The parameter x = A; + 1/2 corresponds to the exponent of the singularity of the above 
density at x = 0, and when x ^ {1/2,3/2,5/2} (resp. x ^ {1;2}) the characteristic 
exponent is given in terms of trigonometric (resp. digamma) functions. 

As a final example we mention general hypergeometric Levy processes, these were 
introduced in [T^ as an example of how to use Vigon's theory of philantropy (cf. [2S]) for 
constructing new Levy processes. They have a characteristic exponent given by 

2 

^(z) = yz^ + iaz + <^i{-iz)<^2{iz), 

where a, a G M and $» is the Laplace exponent of a (possibly killed) subordinator from 
the above mentioned /3-class, which hence takes the form 

$(z) = K + 5z+ ^ (5(1 - a + 7, -7) - 5(1 - a + 7 + z/(3, -7)) 

for a killing rate k > 0. 

These examples illustrate how rich the class of meromorphic Levy processes is. 

Let us now return to working out our algorithm in the case of the American put. 



We start by noting that for the payoff function (4.1) it is straightforward to derive the 
following simplification of the recursion for the functions (V'^ )fc>o stated in Theorem 3.1 

«• 

Proposition 4.3. Suppose that P(Xi < 0) > 0, i.e. X is not a subordinator. Let n > 1. 
There exists a decreasing sequence of points (x^ )fc>o, with x^ := logK, such that for 
each k > 1 we have: 

(i) the point x^f is the unique solution on (— oo,a;^ -^) of the equation in z 



-'^/"E, 



K"), 



Ci(^«(")) -/W = o 
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Kn) 



(ii) V^ is given by 



r(n) 



r(n) 



vir\x) 



J^) 







{^n)V^>y^^on{xl\^). 



Proof. The proof is by induction. First we consider k = 1. Denote for s > 



FAs) 



-r/n^ 



K 



5e^«(n) 



-{K 



(note that this is just the left hand side of the equation in (i) after the substitution 
s = e^). Then Fi(0+) = e'^'/^'K - K < and Fi(exp(4"^)) = Fi{K) > as X not being 
a subordinator imphes that F{X^(„) < 0) > 0. Furthermore, as s H- {K — sa)~^ is convex 
for any a the first term in the formula of Fi is also convex and thus Fi has a unique zero 



on {0,K]. Combining these observations with Theorem 3.1 (i) yields (i) and (ii) above. 
This also shows that V^ (x) > {K — e^) on (x^ , Xq ), while for x G [xq , oo) we have 
Vi {x) > = {K — e^)^ as X not being a subordinator implies ¥x{X^(n) < logi^) > 0, 
so (iii) is also established. 

Now suppose that (i)-(iii) above hold for A; — 1 (IH (i), IH (ii) and IH (iii)). Analogue 
to above, denote for s > the left hand side of the equation in (i) after the substitution 



Fk{s) 



-r/n^ 



v;(!l(iogs + x^(„)) -/(logs). 



By IH (ii) we again have Ffc(0+) = e-'^/'^K - K < 0. Also 



fix 



k-lJ 



-r/n-^_ 



(n) 



yin) 



k~2 



(x 



§(")/ 



< e 



'r/n^_^ 



n) 



r{n) 



Vfc_l(Xg(. 



the equality by IH (i) and the inequality by IH (iii) together with X not being a subordi- 
nator. This translates to Ffc(exp(a;^_]^)) > 0. As above, the first term in the formula of F^ 
is convex since s i— ;■ V"^_{(logs + Q;) is convex for any a, as is clear from the definition (2.4). 



;^(") 



Consequently F^. has a unique zero on (0,K], located in {—oo,x^\), and in combination 



with Theorem 3.ini), statement (i) and (ii) follow. 



It remains to show (iii). It follows from the above and IH (ii) that V,: (x) > V, 



in), 



rin), 



(K-e 
subordinator that 



;^H ;^W 



on (x^',x^_ 



;^(") 



fe-1 



x] 



^]. For X > xj^-^ we have by (ii), IH (ii), IH (iii) and X not being a 



K 



(n), 



X) 



and we are done. 



-^/"E. 



r{n) 



vr-iix.i^^ 



> e 



-'^/"E. 



rin) 



vt'^ix^i.^ 



Vtlix) 



n 
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Note that the above result shows that in the approximating value functions it is 
optimal to stop as soon as X falls below a certain level. Or, equivalently, as soon as the 
stock price S = exp(X) falls below a certain level. This is structurally consistent with the 
exact optimal exercise strategy in the original problem, which is to stop as soon as S falls 
below a certain time dependent boundary (cf. e.g. [12] )■ 



If X satisfies (4.3), i.e. in particular if X is a meromorphic Levy process, then it 



turns out the expectations in the above Proposition |4.3| can be worked out explicitly. The 
proof uses induction over k and consists of tedious yet straightforward computations, it 
is therefore omitted. 



Proposition 4.4. Suppose that X satisfies {4-3) and is not compound Poisson. Fix some 



n 



> 1. For any k > the function V^"^ can be expressed in piecewise form as follows: 



i=o 



-d"'o> 



k-l 



/ Af- k-m 

^.-cro>^.«,,,.,„,.).-.....,„,.)e^.c.».(.„,.)^ 



m=l 



j=0 i=0 

where Xq = logK. Expressions for x^" and all the coefficients involved can be found in 
Appendix\A\ 

Remark 4.5. (i) Note that in the Brownian motion case, the structure of the formulae 
for the V^ 's from Proposition 4-4 matches the corresponding formulae found by 
Carr with his Canadisation method (cf. Section^. That is, the formulae coincide 
except for the values of the coefficients A, B and C , and for the values of the 
boundary points x. 

(a) Any payoff function f that can (piecewise) be expressed as a linear combination of 
functions of the form Ax^e^^ + C for A,B,C G M and z G N can be dealt with in 



similar fashion and yields a similar structure as (4-4)- U o- result as in Proposition 



\4.3\ is not available, i.e. if it is a priori not clear how the maximum in the recursive 



relation (3.1) works out, a way forward is to implement a dynamic decision rule 



where first a formula for 



FJx) 



-^/"E, 



rH 



^fc-l(^5(")) 



is computed (which is still of the same structure as in (4-4V (^"f^d then Vk = max{Ffc, /} 
can be determined in a straightforward way. 
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Finally, a payoff' function that is not of this convenient form mentioned above can 
easily be approximated by a function which is of this convenient form, and it is 
straightforward to estimate the error due to this approximation. 

(Hi) It would be possible to further extend this approach to deal with more complicated 
payoffs involving path dependency, for instance to cases where the payoff at time t 
depends also on the running supremum Xt = sup^<^ Xg as is the case for 'lookback 
options'. A step in the algorithm would then not only require X^(q) but the pair 
{X^(q),X^(q)) rather. Explicit expressions for the law of this pair is available for 
meromorphic Levy processes. 

5 Some numerics for the American put 

In this section we discuss some numerics for the American put, obtained by a computer 



implementation of the result in Proposition 4.4 In particular we present some graphs 



of approximating value functions \4 and some graphs of the approximative optimal 
exercise boundary. For the latter, as is well known (cf. e.g. jl2]) the optimal stopping 



time T* in (1.1) with f{x) = {K — e^)~^ is of the form 



T* = mi{u >0\v{T-u, X„) = /(X„)} = mi{u > | X„ < 6(T - m)} (5.5) 

for some non-increasing function b on ]R>o referred to as the optimal exercise boundary, 
which may hence be defined as 

b{t) := sup{y e R I v{t, y) = f{y)} < log{K) for t > 0. 

Hence, since fo r lar ge n we have V^ ( . ) ~ v{k/n, . ) for all /;; > 1 and recalling the result 
of Proposition 4.3 it should be expectecjj that the set of points {{k/n,x^ )}k>i, in the 



captions of the plots referred to as 'boundary points', form an approximation of the set 
{{k/n,b{k/n))}k>,. 

Next, we briefly describe the two main tasks for creating the graphs of different value 
functions and boundary points presented further below. 

The first step is the calculation of the solutions C+ (0 ^^^ C^ U ) oi z h^ q + '^{iz) 
and of the coefficients c^ (i) and c_ (j) by using Proposition 4.2 for i G {0,...,A^+} 
and j G {0, . . . , A^_}. In case a meromorphic Levy process is used for which iV+ = oo or 
A^_ = oo the infinite sums need to be truncated. For this we used the ad-hoc approach of 
ensuring that the truncated sums represent at least 99% of the mass of X^(„), i.e. (compare 



with (4.3)) we determined M~^,M such that 



^We do not provide a rigorous proof of this, which seems shghtly tricky due to the fact that the 
principle of smooth pasting dictates that at least when the driving Levy process has a Gaussian part v 
'pastes smoothly' onto /, i.e. their first derivatives in the space component coincide. Consequently it is 
not a straightforward application of the Implicit function theorem 



16 



and then normalised the coefficients. It is worth noting that M is larger in the case of 
meromorphic Levy process with finite variation than in the cases with infinite variation. 
Particularly, if cr > the largest part of the probability mass is concentrated on C+(0) 
and C-(0) whereas the other zeroes contribute only a small part. 

The second main step is the recursive computation of the value function and the 



boundary points using formula (4.4). A program for this purpose was written in the 
language C++. We found that in certain cases the algorithm may need more than the 
default 'double precision' as working precision in order to have enough significant digits 
left in the end result. For this we made use of the 'GMP/MPFR/MPFRC++' packages. 
Furthermore we made use of 'OpenMP' to parallelise the operations. 

The Levy processes used in this section are Brownian motion with drift, Brownian 
motion with drift plus exponentially distributed jumps (Kou model) and the /3-class, see 
Section HI 

The Figures [T] and [2] contain different value functions and its boundary points for 
processes of the /3-class. We compare different values of Ai = A2 which are responsible for 
the small jump behaviour. One process is similar to a CGMY process (Ai = 1.8, finite 
variation and infinite activity, see [10]). The other processes have finite variation and finite 
activity (Ai = 0.95) respectively infinite variation and infinite activity (Ai = 2.1). 
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Figure 1: Plots of the approximating value function V^ where the driving process X is of 
the /3-class with parameter values a = 0, ai = 10, 0:2 = 100, /3i = /32 = 0.5, ci = C2 = 100 
and Ai = A2 = 0.95,1.8,2.1. The drift /i is chosen as the solution of ^(— /) = — r. 
Furthermore n = 100, A; = 20, r = 0.05, K = 10000 and M+ = 31, M" = 61 
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X,=V2-1 



0.06 0.08 0.1 0.12 0.14 0.16 0.18 

Time to maturity (T=1) 



Figure 2: Three sets of boundary points obtained using the same parameter values and 
settings as in Figure [T] 

In Figures |3] and |4] the value functions and boundary points for a process of the /3-class 
with infinite variation and infinite activity are plotted for different values of the Gaussian 
coefficient a. If we compare over certain intervals {x^\x^ )^ i > j the curves of the value 
functions and boundary points for different ai and 0"2, cti, (72 G {0, 1, 2, 6} we observe the 
following. If the curve of the value function for cti is steeper on (x^ ,Xj ) than the one 
for a different a2, then the curve of the boundary points for ai appears to be more flat 
on that interval and vice versa. This is in accordance to general optimal stopping theory. 
Note that for a = there seems to be a jump in the boundary at 0, that is Figure |4] 
indicates that 6(0+) < logK if cr = while 6(0+) = logi^ if a > 0. This is consistent 
with the observation in jT^ . 
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10.5 11 



(n) 

Figure 3: Plots of the approximating value function V^ where the driving process X 
is of the /3-class with infinite variation and infinite activity. Parameter values are a = 
0, 1,2,6, ai = 56, 0:2 = 56.4, /3i = (32 = 2,ci = C2 = 1.5 and Ai = A2 = 2.8. Furthermore 
n = 500, k = 100, r = 0.05, K = 10000 and M" = M+ = 1 
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Figure 4: Four sets of boundary points obtained using the same parameter values and 
settings as in Figure |3] 

Finally, Figure [5] illustrates the convergence. Hereby we compare the relative difference 
between the value functions for several n and k with a constant ratio k/n = 0.2. 
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Figure 5: Plots of the relative difference between three approximating value functions, for 
n = 800 and k = 160; n = 200 and k = 40; n = 50 and k = 10. The plots are 'zoomed in' 
to the part of the x-axis where the difference is largest. The case n = 800 and k = 160 
is taken as the base. The underlying /3-class process with infinite variation and infinite 
activity has the following parameters: a = 1,q;i = 1, 0:2 = l,/3i = (^2 = 80, Ci = C2 = 1 
and Ai = A2 = 2.3. The drift fi is chosen as the solution of \E'(— /) = — r. Furthermore 
r = 0.05, K = 10000 and M+ = M" = 1 
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A Formulae for the coefficients appearing in Propo- 



sition 4.4 



This Appendix contains expressions for the coefficients that appear in Proposition 4.4 
recursively in k for some n > 1 fixed. Recall that the starting point for the recursion is 
(4.4) for k = 0, which contains only one coefficient, namely Xq := logK. Throughout 



the sequel, fix some k > 1. 

Let the intervals /^ = [x; , xj^^) for / = 1, . . . , A; — 1, Jq = [xq , oo) and Ik = 

{—oo,Xf^\). Furthermore let F+(x;/,j, /c) be an antiderivative of x i— )■ V^_{(x)e^+ ^^^^ 

on the interval /^ and F_{x;l,j,k) an antiderivative of x ^-)■ V^_^{x)e''- ^^'^ on the 



same interval l\ . Recalling the formula of V^_^ from (4.4), note that -F+(. ;l,j,k) and 
F_{.;l,j,k) consist of sums of antiderivatives of functions of the form x ^-^• x*e^^ for 
i G N, it is hence straightforward to apply partial integration and express them in terms 
of elementary functions. For compactness we do not spell out these formulae in all detail 
here. 



Furthermore, define 



c\'{q) -^-^ &_'{q) 



N+ 



^0 - 2^AnhW^~ 2^ 



in)f\ I T A^ An), 



(h+i) 



T^r\h,j) = (-i)M J2 -?(^)(d"Hg)-cf(j)) 

V 9=0, q^j 

q=0 
( N+ 



q=0 
1=0, QT^J 



<h+l) 



where h eN, j e {0, . . . , N+} for tJ"^ and j GJO, . . . , A^_} for T^"^. 

The coefficients appearing in the formula (4.4) for V^ can now be expressed as follows 



- note that we as usual understand an empty sum as equal to zero, and that the expressions 
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for A^^'{ . , j, . , . ) (resp. Ai" ^( . , j, . , . )) below are valid for all j E {0, . . . , A^+} (resp. 

je{0,...,iV_}): 

i) On the interval [Xk , Xk_i)- 

A^:-\0,j,k,k) = e-^/-c^^\j)[F.{x^-},;k,j,k)-F.{x^-^-0,j,k) 

k-l 

+ 5^(f_(xS;/,j,A;)-F_(x;");/,j,A;))); 

B^''\k,k) = e-'/^To; 
C^''\k,k) = e-"/"K, 

where x^ is the unique solution on {—oo,xj^\) of the equation in z: 

J2 e-^-"'(^)"AL"^(0, q, k, k) - B^''\k, k)e' + C(")(A;, k) = K - e\ 
ii) On the interval [xm ,Xm-i) for m = 1, . . . , /c — 1: 
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k-l 

+ J2 {F^{xP,;l,J,k)~F^{x\-^;l,j,k) 

l=m+l 
k—m~l 

+ e-/n J2 Ai^\h,j,m,k-l)T^r\h,j)h\; 



h=0 
k ~m~ i /'A I • 1 ^ 



A^:\t,j,m,k) = e-'-/" J2 A^:\^ + h~hj,m,k-l)T^;'\h-l,j)- 



h=l 



+ e-"/"^^^AV"H^-l,J,^,'^-l), ViG {l,...,A;-m}; 

V 

A^^\0,j,m,k) = e-^/-c^^\j)(F_{x^:l,;m,j,k)-F_{x^-^;0,j,k) 



m—l 



+ J2{F-i4i;l,J,k)-F4x\''^;l,j,k)] 
1=1 

fc— m— 1 

+ e-/« Y. A^^\Kj,m,k-l)T^^\h,j)h\- 

h=0 
k -rn- i ( h 4- ' — 1^1 



h=l 



k—rn—i 

A^:^\z-l + h,j,m,k-l)T^^\h-l,j)^-^ 
- e-''/"— ^^AL"^(i - 1, j,m, k-l), Vi G {1, . . . , A; - m}; 

C(")(m,A;) = e-"/"C(")(m,A;-l). 
iii) On the interval [xq , oo]: 
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A^:\(},jAk) = e-^/vV")(j)fF^.(4"_\;A;,j,A:)-F+(4");0,j,A;) 

fe-i 

1=1 
fe-i 

+ e-'-/«5^4")(/i-l,j,0,A;-l)TS")(/.-l,j)(/i-l)!; 



^! 

^g-r/„^:+_UZ^(n)^^ - 1, J, 0, A; - 1), V? G {1, . . . , A; - 1}. 
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